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Mean  values  of  the  onshore  and  alongshore  components  of  flow  reflected 
the  influence  of  a rip  current  frequently  observed  just  south  of  the 
instrument  locations.  Strong  harmonics  in  the  spectra  of  sea  surface 
fluctuations  and  particle  velocities  infer  nonlinear  conditions.  Co- 
herence values  between  waves  and  onshore  flow  were  high,  ranging  above 
0.9.  The  coherence  between  waves  and  onshore  flow  was  used  to  separate  th 
turbulence  and  wave-induced  velocity  components.  Over  the  range  of  col- 
lapsing to  spilling  breakers  a reasonable  value  for  the  ratio  of  turbulent 
to  wave-induced  velocity  was  determined  to  be  approximately  0.75.  Satura- 
tion regions  were  found  in  the  wave  and  velocity  energy-density  spectra  at 
higher  frequencies  as  evidenced  by  -5  and  -3  slopes,  respectively.  — - 
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ABSTRACT 


Simultaneous  measurements  of  sea  surface  elevation  and  onshore  and 
alongshore  water  particle  velocities  were  measured  at  three  locations 
within  the  surf  zone  using  two  capacitance  type  pentrating  wave  staffs 
and  three  two-component  electromagnetic  flow  meters.  The  probability 
density  functions,  pdf,  for  the  sea  surface  elevation  were  always  highly 
positively  skewed,  whereas  the  pdf's  for  the  velocities  were  both  nega- 
tively and  positively  skewed.  Mean  values  of  the  onshore  and  alongshore 
components  of  flow  reflected  the  influence  of  a rip  current  frequently 
observed  just  south  of  the  instrument  locations.  Strong  harmonics  in 
the  spectra  of  sea  surface  fluctuations  and  particle  velocities  infer 
nonlinear  conditions.  Coherence  values  between  waves  and  onshore  flow 
were  high,  ranging  above  0.9.  The  coherence  between  waves  and  onshore 
flow  was  used  to  separate  the  turbulence  and  wave-induced  velocity 
components.  Over  the  range  of  collapsing  to  spilling  breakers  a rea- 
sonable value  for  the  ratio  of  turbulent  to  wave-induced  velocity  was 
determined  to  be  approximately  0.75.  Saturation  regions  were  found  in 
the  wave  and  velocity  energy-density  spectra  at  higher  frequencies  as 
evidenced  by  -5  and  -3  slopes,  respectively. 
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I.  INTRODUCTION 

A.  HISTORICAL  PERSPECTIVE 

The  surf  zone  is  an  area  bounded  on  the  seaward  side  by  the  point 
where  waves  first  begin  to  break,  and  on  the  landward  side  by  the  point 
of  maximum  run-up  on  the  beach  slope.  Although  beach  profiles  and  waves 
within  the  surf  zone  have  been  studied  by  many,  the  kinematics  of  wave 
forms  and  water  particle  velocities  in  the  zone  have  remained  something 
of  a mystery. 

The  problems  encountered  have  been  practical  as  well  as  theoretical. 
Wave  theories  developed  for  deep  water  Waves  do  not  correctly  characterize 
the  motion  which  occurs  as  the  waves  break.  Direct  measurements  of  break- 
ing waves  appear  to  be  the  most  viable  means  of  approaching  the  problem. 
However,  this  has  been  hampered  by  the  surf  zone's  hostile  environment, 
which  is  difficult  to  reproduce  in  the  laboratory,  and  by  inadequate 
instrument  design.  Improvements  in  the  latter  area  has  resulted  in  sturdy, 
sensitive,  but  expensive  measuring  devices  which  have  rapid  response  times. 

The  earliest  study  of  significant  importance  was  conducted  by 
Iversen  (1953)  in  which  he  used  photographic  techniques  to  obtain  a 
Lagrangian  description  of  water  particle  motion  under  breaking  waves. 

The  laboratory  channel  limited  the  wave  type  to  plunging  and  surging 
breakers. 

Inman  (1956)  was  one  of  the  first  to  make  sophisticated  field 
measurements,  measuring  the  drag  force  on  a cantilevered  sphere  to  infer 
water  particle  motion.  Miller  and  Zeigler  (1964)  measured  in  situ 
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particle  motions  using  acoustic  and  electromagnetic  current  meters, 
then  compared  their  findings  with  higher  order  wave  theory  and  found 
some  qualitative  agreement.  Walker  (1969)  made  similar  measurements 
using  propeller  type  flow  meters.  Meadows  (1976)  used  ducted  impeller 
flow  meters  at  equally  spaced  vertical  positions  in  order  to  measure 
longshore  components  of  flow.  Huntley  (1976)  utilized  a single  two- 
} component  electromagnetic  flow  meter  in  an  effort  to  obtain  a value 

of  the  friction  coefficient.  Bowen  and  Huntley  (1974)  made  measure- 
ments of  the  nearshore  velocity  fields  using  as  many  as  three  two- 
component  electromagnetic  flow  meters.  FLihrboter  and  Biisching  (1974) 

: measured  simultaneous  orbital  velocities  and  water  levels  using  a two- 

^ • component  current  meter  and  two  pressure  type  wave  meters.  Thornton 

• (1969),  Steer  (1972),  Thornton  and  Richardson  (1973),  Bub  (1974)  and 

* Galvin  (1975)  used  pressure  meters,  capacitace  wave  gauges  and  electro- 

I magnetic  current  meters  to  compute  particle  velocities  and  surface  profiles 

i 

I within  the  surf  zone.  The  work  presented  herein  is  an  extension  of  these 

latter  studies. 

B.  OBJECTIVE 

The  objective  of  this  research  is  to  study  the  kinematics  of  water 
\ particles  within  breaking  waves  and  within  the  surf  zone.  In  particular, 

I the  experiment  was  performed  on  a flat  beach  in  order  to  measure  the 

characteristics  of  spilling  breakers.  Simultaneous  measurements  were 
. made  of  the  instaneous  sea  surface  elevation  and  two  orthogonal  water 

particle  velocities,  onshore  and  alongshore,  at  three  fixed  locations 
in  the  surf  zone.  Estimates  of  the  probability  density  functions  and 


power  spectra  of  the  wave  heights  and  particle  velocities  are  made.  The 
suitability  of  using  linear  wave  theory  as  a spectral  transfer  function 
in  determining  velocity  spectral  components  from  the  power  spectrum  of 
the  waves  is  measured.  Computed  velocity  spectra  are  compared  with 
actually  measured  velocity  spectra.  Coherence  between  waves  and  onshore 
flow  is  used  to  separate  the  turbulence  and  wave-induced  velocity  com- 
ponents in  order  to  determine  the  ratio  of  turbulent  to  wave-induced 
velocity.  Spectral  estimates  are  analyzed  to  ascertain  the  slope  of 
the  saturation  region  at  higher  frequencies. 


II.  MEASUREMENTS 


A.  EXPERIMENT 

The  experimental  site  is  just  north  of  La  Jolla,  California.  Data 
presented  in  this  study  was  taken  as  part  of  a larger  scale  experiment 
conducted  from  15  February  to30  March  1977  in  the  Southern  California 
bight  area.  The  over-all  purpose  of  the  experiment  was  to  test  and 
evaluate  the  instrumentation  package  which  was  selected  for  use  aboard 
SEASAT.  The  apparatus,  which  will  measure  waves,  winds,  and  sea  surface 
temperature,  was  being  flown  on  board  NASA  airplanes.  A major  ground 
truth  program  was  concurrently  being  conducted  at  Torrey  Pines  Beach. 

Measurements  of  waves,  currents,  and  sediment  transport  ac  Torrey 
Pines  Beach  were  made  using  the  following  equipment.  A five  element 
linear  array  of  pressure  sensors  was  positioned  offshore  along  the  ten 
meter  contour  to  measure  wave  energy  and  direction.  Located  in  the  surf 
zone  inshore  of  the  array  were  three  support  towers  in  a line  perpendi- 
cular to  the  beach.  A surface  piercing  capacitance  wave  gauge  and  a two- 
component  electromagnetic  flow  meter  were  mounted  on  each  tower.  A 
nephelometer  was  installed  on  the  middle  support  tower.  To  the  north 
of  the  three  towers  was  an  array  of  four  flow  meters  which  was  used  to 
delineate  the  areal  and  temporal  variation  in  alongshore  currents. 

Figure  1 shows  instrument  locations.  Resistance-wire  run-up  meters  were 
installed  in  the  swash  zone  at  four  locations.  All  instrumentation 
was  battery  operated. 

The  offshore  pressure  sensor  outputs  were  telemetered  directly  to 
the  shore  processes  lab  at  Scripps  Institute  of  Oceanography,  which  is 
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located  one  mile  to  the  south.  All  data  from  the  other  instruments 
were  cabled  to  one  of  two  transmitting  terminals  on  the  beach  where  it  was 
then  telemetered  to  the  shore  processes  lab  and  recorded.  All  the  equip- 
ment except  the  nephelometer  were  in  position  prior  to  high  tide  on  9 
March.  The  nepholometer  was  installed  on  14  and  15  March.  Data  was 
accumulated  during  high  tide  on  9,  10,  11,  16,  17,  18,  19,  21  and  23 
March. 

Several  other  smaller  scale  experiments  were  conducted  at  the  same 
time.  At  an  elevation  of  approximately  300  feet  on  a cliff  overlooking 
Torrey  Pines  Beach  a radar  installation  was  established.  The  radar 
images  of  waves  were  then  photographed  and  analyzed. 

Sediment  transport  experiments  were  conducted  on  11,  21,  and  23 
March.  Bed  load  transport  was  measured  using  the  methods  described  by 
Inman  end  Komar  (1970)  to  trace  the  movement  of  fluorescent  dyed  sand. 
Suspended  sediments  were  measured  in  situ  by  swimmers  using  a mechanical 
water  sampling  devide.  These  samples  were  taken  along  with  nepholometer 
readings  in  order  to  compare  methods.  During  these  sediment  transport 
investigations,  Lagrangian  floats  made  from  wine  bottles  weighted  by 
sand  were  used  to  determine  average  long-shore  current  speeds. 

B.  INSTRUMENTATION. 

The  manner  in  which  waves  break  depends  very  much  on  the  ;haracter- 
istics  of  the  beach  and  near  shore  bottom  slope.  (Table  I)  Spilling  breakers 
occurred  most  frequently  at  the  Torrey  Pines  Beach  site.  Figure  2 shows 
a typical  beach  profile  and  the  location  of  the  instruments  which  provided 
the  data  for  this  paper. 
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TABLE  I.  BEACH  AND  WAVE  CHARACTERISTICS 

LOCATION 
DATE 

BEACH  SLOPE 
SAND  TYPE 

SAND  MEDIAN  GRAIN  SIZE 
BREAKER  TYPE 
WAVE  PERIOD 
WAVE  HEIGHT 


TORREY  PINES  BEACH 
9-23  MARCH  1977 
FIGURE  2 and  APPENDIX  A 
QUARTZ 

I. 2  MILLIMETERS 
SPILLING 

II. 9  to  15.9  SECONDS 
UP  TO  1.5  METERS 
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1 . Wave  Gauge 


The  wave  gauges  were  of  the  capacitance  type  fashioned  from 
3/8  inch  outside  diameter  stainless  steel  rod.  The  rod  was  tightly 
covered  with  1/16  inch  wall  thickness  polypropylene  tubing.  This  li- 
near, highly  sensitive  instrument  has  proven  to  be  sufficiently  sturdy 
to  withstand  the  rigors  of  the  unfriendly  surf  zone.  The  gauges  operate 
on  the  principle  that  a change  in  the  plate  dimension  of  the  capacitor 
changes  its  capacitance  and  consequently  the  circuitry  output  voltage. 

In  these  gauges,  the  insulated  steel  rod  and  sea  water  act  as  the  plates 
and  the  insulation  functions  as  the  dielectric.  As  the  water  level  flue 
tuated,  the  capacitance  changed.  Fluctuations  were  sensed  by  a transis- 
torized circuit  powered  from  the  beach.  The  circuit  was  designed  by 
McGoldrick  (1969).  The  electronics  packaged  was  housed  in  a watertight 
brass  case  which  was  mounted  on  the  tower.  This  allowed  the  connecting 
leads  to  be  relatively  short,  hence  minimizing  wire-to-wire  capacitance. 
All  gauges  were  statically  calibrated  in  the  laboratory  prior  to  the  ex- 
periment. Accuracy  was  estimated  to  be  -.005m.  The  calibration  plots 
are  shown  in  Appendix  B. 
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2.  Flow  Meter 

The  flow  meters  were  Marsh-McBirney  model  511  electromagnetic 
water  current  meters.  Operation  of  the  meter  is  based  on  Faraday's  prin- 
ciple of  electromagnetic  induction.  The  meter  measures  water  particle 
velocity  in  two  orthogonal  directions  through  a range  of  zero  to  three 
m/sec  with  a maximum  output  error  of  two  percent  of  full  scale  reading. 
The  sensor  has  a variable  time  constant,  a setting  of  0.2  sec.  was  used. 

The  flow  meters  were  dynamically  calibrated  using  the  method  of 
Thornton  and  Krapohl  (1974).  This  procedure  involves  utilizing  an  os- 
cillating platform  driven  by  a variable  speed  motor  geared  to  an  eccen- 
tric throw  arm.  This  technique  was  employed  in  order  to  determine  the 
meter's  characteristics.  Measurement  accuracy  was  determined  to  be 
*.005  m/sec.  The  only  difficulty  encountered  in  their  utilization  oc- 
curred when  Flow  Meter  #1  was  buried  on  19  March  due  to  the  heavy  sedi- 
ment transport. 

The  measuring  instruments  were  attached  to  a 3.6m  high  tower  con- 
structed of  steel  pipe  with  an  outside  diameter  of  6. 3cm. A 0.5m  diameter 
base  plate  was  placed  about  0.6m  from  the  bottom  of  the  tower.  This 
configuration  allowed  the  tower  to  be  sunk  approximately  0.7m  into  the 
sand.  The  towers  were  also  supported  by  steel  guy  wires  fastened  to 
blade  anchors  driven  into  the  sand.  Thus, tower  movement  and  vibration 
were  negligible. 

The  towers  were  placed  on  a line  perpendicular  to  the  shore  at 
specified  distances,  and  were  erected  during  low  tide  when  the  beach 
was  more  accessible.  As  indicated  earlier,  measurements  were  taken  at 
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high  tide.  Instruments  were  arranged  so  as  to  be  in  the  same  vertical 
plane.  A carpenter's  level  was  used  to  establish  axis  alignment  with  an 
estimated  error  of  *2  degrees.  A typical  tower  and  instrument  arrangement 
is  shown  in  Figure  3. 
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Schematic  of  Tvpical  Instrument  Tower  with  Sensors 


I 

III.  ANALYSIS  OF  DATA 

The  data  were  digitized  at  a rate  of  four  samples  per  second, 
corresponding  to  a sample  interval  of  0.25  seconds.  This  resulted  in  a 
Nyquist  frequence  of  2.0  Hz.  The  sampling  rate  was  considered  suffi- 
ciently high  enough  to  avoid  aliasing  of  energy  into  the  portion  of  the 
spectra  which  was  of  interest. 

Continuous  time  series  records  of  all  measurements  were  taken  for 
roughly  2-1/2  hours  each  day.  Record  lengths  of  24  minutes  from  each 
data  set  were  subjected  to  analysis.  Criteria  which  determined  record 
lengths  included  economy  of  computer  usage  and  resolution  over  the  fre- 
quency range  of  interest. 

A maximum  lag  time  was  chosen  as  5%  of  the  record  length.  This  re- 
sulted in  a spectral  band  width  resolution  of  0.007  Hz  and  each  spectral 
estimate  having  40  degrees  of  freedom.  The  ninety  percent  confidence 
limits  for  40  degrees  of  freedom  using  a chi-square  distribution  are  found 
to  be  between  0.72  and  1.51  of  the  measured  power  spectral  estimates. 

A mean  value  was  computed  for  all  data  sets  and  the  data  was  linearly 
detrended  to  remove  tidal  effects.  Variance,  standard  deviation  and 
average  period  were  calculated.  The  average  period  was  determined  by 
calculating  the  time  between  zero  upcrossings.  The  calculated  average 
periods  are  lower  than  the  visually  observed  periods  due  to  a number  of 
perturbations,  such  as  noise  and  capillary  waves,  increasing  the  zero 
upcrossing  occurrences.  A probability  density  function  for  each  data  set 
was  calculated  and  plotted.  Comparisons  were  made  with  Gaussian  and  Gram- 
Charlier  distributions  using  the  chi-square  goodness-of-fit  test.  Variance, 
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standard  deviation,  skewness  and  kurtosis  of  the  distributions  were 
computed  and  are  summarized  in  Table  II. 

For  each  data  set  an  auto-covariance  function  was  determined  and  then 
smoothed  with  a Parzen  window.  Applying  a Fourier  transform  to  the  smoothed 
auto-covariance  function  resulted  in  the  power  spectrum.  An  examination  of 
the  power  spectrum  shows  the  regions  of  greatest  potential  or  kinetic 
energies  and  their  respective  frequencies. 

Cross-spectra  were  computed  by  calculating  a cross-covariance  between 
data  sets,  smoothing  it  with  a Parzen  window  and  applying  a Fourier  trans- 
form. Coherence  and  phase  were  then  determined  from  the  cross-spectrum. 

A coherence  and  phase  versus  frequency  plot  will  indicate  the  regions 
and  degree  of  linear  relationship  and  phase  between  two  data  sets. 


TABLE  II.  COMPUTED  STATISTICS 
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TABLE  II.  COMPUTED  STATISTICS  (Continued) 

STANDARD  CHI-SQUARE  PARAMETER 

DATE  INSTRUMENT  MEAN  VARIANCE  DEVIATION  SKEWNESS  KURTOSIS  GAUSSIAN  GRAM-CHARLIER 
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Mean  and  variance  for  the  wave  heights  is  given  in  meters  squared;  mean  and  variance  for 
the  velocity  components  is  given  in  meters  squared  per  second  squared. 

Standard  deviation  for  the  wave  heights  is  given  in  meters;  standard  deviation  for  the 
velocity  components  is  given  in  meters  per  second. 


IV.  RESULTS 

A.  QUALITATIVE  DESCRIPTION 

Observation  of  various  breaker  types  which  occur  on  different  beaches 
has  resulted  in  a number  of  similarities  being  noted.  A typical  analog 
record  of  waves  and  onshore  velocities  obtained  at  Torrey  Pines  Beach  is 
shown  in  Figure  4.  In  general,  a sawtooth  shaped  profile  is  created  when  ' 

there  is  a quick  drawdown  of  water  just  before  the  breaker  arrives,  fol-  I 

j 

lowed  by  a steep,  vertical  leading  edge,  and  a sloping  profile  toward  | 

the  trailing  edge.  On  the  trailing  edge,  secondary  waves  are  often  ob-  | 

served.  Secondary  waves  are  harmonics  of  the  primary  wave  frequency  and 

denote  strongly  nonlinear  waves.  At  Torrey  Pines  Beach  the  breaking  ; 

waves  are  generally  of  the  spilling  variety  which  spill  rapidly  at  the 

crest  and  move  down  the  wave.  j 

I 

B.  MEAN  VALUES 

A mean  value  was  computed  for  all  data  sets  and  linearly  detrended 
to  remove  tidal  effects.  The  sea  surface  elevation  mean  values  ranged 
from  0.67  to  1.33.  Wave  heights  increased  till  the  18th  and  then  de- 
creased during  the  remaining  period  data  was  taken,  reflecting  well  the 
actual  environment. 

Analysis  of  the  mean  values  for  both  the  onshore  and  alongshore 
components  of  flow  was  complicated  by  an  apparent  rip  current.  A rip 
current  was  frequently  observed  located  just  south  of  the  array  of  wave 
gauges  and  flow  meters.  The  on-offshore  flow  for  the  five  days  under 
consideration  at  Flow  Meter  #1  was  always  directed  offshore  at  about  0.1 
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ALONGSHORE  VELOCITY 


Figure  4.  Typical  Analog  Record  of  Waves  and 
Velocities  from  Torrey  Pines  Beach. 


m/sec.  At  Flow  Meter  #2  the  horizontal  water  particle  velocity  was  di- 
rected onshore  at  roughly  0.25  m/sec.  The  flow  at  Flow  Meter  #3  was 
away  from  the  coast  on  16,  17,  and  18  March  at  slow  speeds  and  towards 
the  shore  on  21  and  23  March  moving  very  slowly.  Concurrently,  the  long- 
shore flow  at  Flow  Meters  #1,  #2  and  #3  was  generally  moving  downcoast 
with  the  peak  velocity  varying  from  flow  meter  to  flow  meter  from  day  to 
day.  Similar  results  were  reported  by  Huntley  and  Bowen  (1974)  for  a 
nearshore  circulation  cell  which  they  attributed  to  a rip  current. 

C.  PROBABILITY  DENSITY  FUNCTIONS 

1 . • Gaussian  and  Gram-Charl ier  Frequency  Distributions 

The  computed  probability  density  function  was  compared  to  the 
Gaussian  and  Gram-Charl ier  distributions  and  tested  for  comparison  using 
the  chi-square  goodness-of-fit  test.  The  closer  the  chi-square  fit 
parameter  is  to  zero,  the  better  the  fit  (Table  II). 

The  Gaussian  probability  density  function  is  given  by 

PGa**’  = — ! 

a yTn 

where  a is  the  standard  deviation  and  p is  the  mean.  The  Gaussian 
distribution  is  completely  described  by  the  mean  and  variance,  and  has 
zero  skewness  and  a kurtosis  equal  to  three.  This  approximates  the 
sea  surface  in  deep  water  conditions.  However,  nonlinearities  found  in 
the  surf  zone  introduce  skewness  and  kurtosis  values  that  deviate  from 
Gaussian  and  result  in  a distribution  more  closely  approximated  by  the 
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Gram-Charl ier  distribution.  This  is  to  be  expected  in  that  the  Gram- 
Charlier  pdf  is  calculated  using  the  assigned  parameters  from  the  data, 
and  includes  moments  higher  than  the  second.  The  frequency  distribution 
of  Wave  Gauge  #1  on  23  March  is  shown  in  Figure  5.  Appendix  D contains 
distributions  calculated  for  other  data  sets. 

Distributions  for  all  flow  velocities  are  either  unimodal  or  have 
less  pronounced  secondary  peaks  compared  to  the  distributions  of  wave 
heights.  The  distribution  of  the  onshore  flow  of  Flow  Meter  #1  on  23 
March  is  shown  in  Figure  6.  As  noted  in  Table  II,  this  results  in  nearly 
equal  fit  parameters  with  the  Gram-Charl ier  distribution  being  generally 
smaller  and  thus  a better  fit.  Fit  parameters  for  velocity  components 
are  much  closer  to  zero  than  are  the  values  calculated  for  wave  components 
since  the  flow  meters  do  not  experience  the  surface  irregularities  to 
the  same  degree  as  do  the  wave  gauges. 

2.  Skewness  and  Kurtosis 

The  sea  surface  distribution  was  found  to  be  positively  skewed. 

This  means  that  positive  values  are  large,  but  less  frequent  than  negative 
values  which  are  more  frequent,  but  smaller,  indicating  a greater  amount 
of  time  below  the  mean  water  level.  In  other  words,  the  crests  are  steeper 
and  more  peaked,  whereas  the  troughs  are  elongated  and  flatter  which  is 
how  waves  in  shallow  water  theoretically  appear. 

It  should  follow  that  the  wave-induced  particle  velocities  would 
have  distribution  similar  to  the  waves.  However,  of  the  34  cases  listed 
in  Table  II,  17  had  positive  skewness  and  17  had  negative  skewness  with  no 
readily  apparent  pattern.  Skewness  of  the  onshore  flow  at  Flow  Meter  #1 


was  always  positive.  Alongshore  flow  at  the  same  location  had  negative 
skewness  in  three  of  four  cases.  Onshore  water  particle  velocity  at 
Flow  Meter  #2  had  negative  skewness  for  all  six  days.  The  skewness  of 
alongshore  flow  at  the  same  point  was  positive  in  five  of  six  cases. 

The  onshore  component  of  flow  at  Flow  Meter  #3  had  a positive  skewness 
for  16,  17,  and  18  March  and  a negative  skewness  for  the  following  three 
days  on  which  data  was  taken.  Skewness  of  the  alongshore  flow  at  Flow 
Meter  #3  was  positive  for  four  of  the  six  days.  The  reason  for  this 
disparity  is  not  known.  Similar  anomalous  results  have  been  found  by 
Thornton  and  Galvin  (1975). 

Kurtosis  indicates  the  peakedness  of  each  parameter.  Just  prior 
to  breaking  the  waves  achieve  the  greatest  degrees  of  peakedness.  Visual 
observations  indicated  that  most  of  the  waves  were  breaking  at  or  near 
Wave  Gauge  #1  on  17,  18,  and  19  March;  this  correlates  with  their  high 
kurtosis  values  on  these  days.  Values  of  kurtosis  for  all  flow  veloci- 
ties again  had  no  apparent  pattern  with  values  ranging  from  2.6  to  8.6. 
The  reason  for  this  is  not  known. 

D.  SPECTRAL  ANALYSIS 

1 . Sea  Surface  Elevation  and  Water  Particle  Velocity 

The  power,  coherence,  and  phase  spectra  of  Wave  Gauge  #1  and  the 
onshore  component  of  the  flow  of  Flow  Meter  #1  on  23  March  are  shown  in 
Figure  7.  The  spectrum  of  the  wave  surface  elevation,  which  is  char- 
acteristic of  all  the  records,  show  a narrow  banded  peak  at  a frequency 
of  0.07  Hz,  corresponding  to  a period  of  14.3  sec.  Also  evident  are 
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peaks  at  0.14  and  0.21  Hz  which  appear  to  be  harmonics  of  the  primary 
peak  at  0.07  Hz.  These  harmonics  appear  to  be  physical  as  was  observed 
from  the  strip  chart  earlier.  The  harmonics  probably  have  some  energy 
contribution  due  to  the  Fourier  computational  technique.  Appendix  E 
contains  spectra  calculated  for  all  other  wave  and  onshore  flow  data  sets. 

The  coherence  values  were  high,  ranging  to  greater  than  0.9  in 
the  maximum  energy  portion  of  the  wave.  Coherence  begins  to  fall  off 
at  about  0.5  Hz  except  on  the  17th  and  18th  of  March,  the  two  days  on 
which  short-crested  waves  were  noted.  On  these  days  the  coherence  begins 
to  decrease  at  about  0.3  Hz.  The  decrease  in  coherence  is  probably  due 
to  the  relative  increase  in  turbulence,  the  general  decrease  in  energy 
level,  and  wave  energy  spreading.  The  decrease  in  coherence  due  to  wave 
energy  spreading  is  discussed  in  Section  E. 

The  phase  angle,  which  according  to  linear  theory  should  be  zero, 
had  an  average  value  of  less  than  20  degrees  over  the  highly  coherent 
band  of  prominent  wave  energy.  The  phase  angle  tends  to  increase  slightly 
over  the  coherent  band  and  becomes  random  for  the  noncoherent  region  of 
the  spectrum. 

In  order  to  compare  waves  and  flow  velocity,  the  wave  profile 
spectrum  was  converted  to  a theoretical  velocity  spectrum  for  comparison 
with  the  measured  flow  velocity  spectrum  using  the  linear  transfer  function 
such  that 

Su(f)  = |H(f)|  ^ (f) 

The  transfer  function  is  given  from  linear  wave  theory 


(2) 


I 

I 


\ 


I 


I 


I 

H(f)  = 2 irf  coshk(h-<-z) (3) 

sinh  kh 


where 

k is  the  wave  number, 

h is  the  mean  water  depth,  and 

z is  the  depth  of  flow  meter  below  the  mean  depth  (m) 

Quantitatively,  the  results  showed  that  linear  theory  overestimated 
wave-induced  velocity  spectral  components  by  about  50%;  in  previous 
studies  by  Thornton  (1976)  where  the  waves  were  long  crested  and 

arrived  at  near  normal  incidence,  the  opposite  was  true,  in  that  linear 
theory  underestimated  the  magnitude  of  the  wave  spectrum.  The  difference 
in  results  is  probably  due  to  wave  directionality,  since  when  the  wave 
approaches  at  any  angle  other  than  normal  the  flow  meter  measures  only  a 
component  of  the  actual  flow. 

2.  Onshore  and  Alongshore  Flow 

The  cross-spectra  between  the  onshore  and  alongshore  components 
of  flow  were  calculated  with  the  intent  of  finding  the  relation  between  these 
two  parameters.  Power,  coherence,  and  phase  spectra  of  the  onshore  and 
alongshore  flows  at  Flow  Meter  #3  on  23  March  are  shown  in  Figure  8. 

The  spectra,  which  are  representative  of  all  the  records,  shows  very 
little  coherence  between  velocity  components.  The  phase  angle  is  ran- 
dom over  the  noncoherent  band. 

The  reason  for  the  lack  of  coherence  between  the  two  orthogonal 
flow  components  is  believed  to  be  due  to  wave  directionality.  As  was 
seen  earlier  the  onshore  component  of  flow  appears  highly  wave-induced; 
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concomitantly,  the  alongshore  flow  component  is  at  near  normal  incidence 
to  the  wave  direction  and  thus  has  only  a small  wave-induced  contribution. 
The  variance  of  the  alongshore  component  of  flow  is  significantly  less  than 
the  on-offshore  component  except  at  zero  frequency,  which  reflects  the 
mean  longshore  current  variation. 

E.  TURBULENT  VELOCITY  VERSUS  WAVE- INDUCED  VELOCITY 

The  total  velocity  can  be  separated  into  components  of  a mean,  plus 
wave-induced,  plus  turbulent  velocity 

u = u + U + u’ . (4) 

The  wave-induced  and  turbulent  velocity  spectral  components  are  assumed 
to  be  statistically  independent.  For  unidirectional  waves  aligned  with 
the  measured  velocity  component,  the  co-spectra  between  waves  and  velocity 
is  then  given  by 

W' ' Wf’- 

Assuming  statistical  independence,  the  horizontal  velocity  spectrum 
is  obtained  from 


Sjf)  = S^.(f)  + Su(f).  (6) 

A further  assumption  is  made  that  the  waves  and  wave-induced  velocities 
are  given  by  a constant  parameter  linear  process  where  the  coherence 
is  equal  to  unity 


’'toU) 


'U?7(f) 


= 1. 


(7) 


Substituting  (5),  (6),  and  (7)  into  the  definition  of  coherence  between 
the  total  horizontal  velocity  and  waves  results  in 
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y (f)  =! 
ury  ' I 


1 + Vm 

V (f) 


-1 


= 


(8) 


Increasing  lack  of  coherence  is  due  to  an  increasing  high  ratio  of  turbu- 
lence. The  coherence  values  indicate  the  percent  of  total  velocity 
which  is  associated  with  the  wave  velocity. 

Wave-induced  velocity  can  be  determined  using  (8) 


Su(f)  (f)s.,(f). 


nr} 


(9) 


The  wave-induced  velocity  spectral  component  calculated  in  this  way 
should  be  less  than  the  actual  value,  the  reason  for  this  being  that 
the  linear  coherence  between  waves  and  velocities  will  always  be  under- 
predicted due  to  nonlinearities  which  are  always  large  in  a breaking 
wave  and  also  as  a result  of  directional  spreading  of  the  incident 
wave  energy  (Battjes  (1974)  ). 

The  turbulent  velocity  spectrum  is  derived  by  subtracting  the  cal- 
culated velocity  spectrum  from  the  measured  velocity  spectrum, 

Su.(f)  = - Su(f).  (10) 

The  ratio  of  the  turbulent  velocity  intensity  to  the  wave-induced 
velocities  intensity  is  shown  in  Figure  9.  The  data  used  for  collapsing 
and  plunging  breakers,  denoted  by  squares,  was  measured  during  earlier 
studies  by  Thornton  et  (1976).  The  collapsing  and  plunging  breakers 
were  measured  on  beaches  where  the  waves  were  long  crested  and  arrived 
at  near  normal  incidence.  The  stars  represent  data  collected  during  this 
experiment.  A value  of  the  velocity  intensity  is  given  by  taking  the 
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turbulent  to  wave- induced 


square  root  of  the  variance  which  is  computed  by  integrating  the 


velocity  spectrum,  across  all  frequency  bands 

The  measurements  show  a small  spread  of  values  for  the  ratio  of  turbulent 
to  wave-induced  velocity  intensities.  The  average  ratio  calculated  is 
approximately  0.82.  This  would  indicate  that,  to  a first  approximation, 
the  wave-induced  kinetic  energy  is  approximately  1.2  times  the  turbulent 
kinetic  energy.  Thus  the  velocities  in  the  surf  zone  on  the  average  are 
highly  wave-induced. 

If  waves  are  not  unidirectional,  and  generally  they  are  not,  it  is 
necessary  to  account  for  directional  spreading.  The  wave  angle  striking 
the  wave  staff  and  the  type  of  wave,  long  or  short  crested,  must  be  con- 
sidered. In  general,  short  crested  waves  predominated  at  Torrey  Pines 
Beach.  On  the  17th  and  18th  of  March  they  were  particularly  noticeable. 

The  coherence  between  surface  elevation  and  velocity  components  of  short 
crested  waves  is  less  than  that  of  long  crested  waves  due  to  directional 
spreading.  Using  the  wave  train  model  developed  by  Yefimov  and  Khristoforov 
(1971)  the  ratio  of  coherence  between  waves  and  onshore  flow  for  the  short 
crested  waves  to  that  for  long  crested  waves  is  given  by 
^ (1  * l/3ctnV  )-' 

( 1 * S„,{f)/Sy(f)[l/cos^  ] )■' 

The  model  describes  long  crested  waves  using  a train  of  two-dimensional, 

planar  random  waves  traveling  in  direction  . This  corresponds  to  a I 

bivariate  spectrum  given  by  I 

i 

i 

y 
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(13) 


4<^  (f,^l)  = S^(f)6(^-^o)  , 

where  5«^-^o  ) is  the  Dirac  delta  function. 

Short  crested  waves  appear  as  a train  of  three-dimensional  waves 

whose  average  direction  of  spectral  components  is  assumed  to  be  of  the  form 

(2/t7)  S^(f)  cos'^ie-Oo)  , 

0 fori  0 - ^ol>  TT  /2  , (14) 

Neglecting  turbulence  in  order  to  examine  the  effect  of  directionality, 

the  ratio  of  the  coherence  between  waves  and  wave-induced  velocities  for 

short  crested  waves  to  that  of  long  crested  waves  becomes 


'Pr.  (I'.^S) 


(1  + l/Sctn^^o)’"' 


(15) 


As  an  example, when  Oo  equals  an  angle  of  30°,  the  ratio  equals  0.9. 


Conversely,  the  influence  of  turbulence  on  the  ratio  of  coherence 
between  waves  and  velocities  for  short  crested  waves  to  that  of  long 
crested  waves  can  be  examined  by  disregarding  the  directionality  and 
letting  ^o=0.  The  ratio  is  given  by 

y^Jf.^s)  = ( 1 + 4s  .(f)/3Su(f)  )■’ 

f . (16) 

y^f.^l)  ( 1 + Sj^.(f)/Su(f)  )"' 

A summary  of  the  ratios  considering  the  turbulence  and/or  direction- 
ality factors  at  various  angles  is  tabulated  in  Table  III.  The  turbulence 
factor,  Sy.(f)/Sy(f ),  is  assumed  for  comparison  purposes  to  be  equal  to 
0.75,  the  average  value  for  the  measured  data  in  Figure  9 as  denoted  by 
the  squares. 


The  coherence  between  waves  and  velocities  was  used  to  calculate  the 

wave-induced  velocity  (equation  11).  This  method  under-predicts  the  actual 

wave-induced  velocity  because  of  nonlinearities  and  wave  energy  spreading. 

Since  the  turbulent  velocity  is  computed  by  subtracting  the  calculated 

wave-induced  velocity  spectrum  from  the  measured  velocity  spectrum 

(equation  12)  the  smaller  wave-induced  velocity  value  makes  the  turbulent 

velocity  value  larger.  The  average  ratio  of  the  turbulent  velocity 

intensity  to  the  wave-induced  velocity  intensity  was  calculated  to  be 

0.825.  In  order  to  correct  for  the  directional  spreading  it  is  necessary 

2 2 

to  multiply  by  the  ratio  V ^^)/y  (f.^l)-  As  an  example,  the 

ratio  of  the  coherences  for  an  angle  of  15°  is  equal  to  0.891.  Multiplying 
this  value  by  the  average  ratio  results  in  a turbulence  to  wave-induced 
ratio  of  0.75,  which  is  equal  to  the  earlier  measured  mean  value.  Hence 


over  the  range  of  collapsing  to  spilling  breakers  a reasonable  value  for 


i 


► 
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approximately  0.75. 

F.  SATURATION  REGION  IN  THE  SPECTRUM  OF  BREAKING  WAVES 

Waves  on  Torrey  Pines  Beach  characteristically  broke  as  spilling 
breakers  with  the  crest  sliding  down  the  face  of  the  wave.  When  the 
fluid  particles  at  the  free  surface  move  forward  at  a greater  speed 
than  the  wave  speed,  c,  breaking  occurs.  Strong  harmonics  observed 
in  the  wave  and  velocity  spectra  (Figure  7)  are  indicative  of  energy 
being  transferred  from  low  to  higher  frequencies.  The  energy  is  even- 
tually dissipated  by  viscosity  at  the  highest  frequencies.  When  the 
transfer  of  energy  is  not  rapic^^enough  to  balance  the  increase  in  energy 
density  of  the  waves  during  s,Hoaling,  breaking  occurs.  The  waves  are 

said  to  be  "saturated"  with'energy.  Hence,  a regiop'  of  saturation 

J f 

would  be  expected  through' which  energy  is  transferred  from  the  low  to 
higher  frequencies.  During  shoaling  and  breaki/g  of  waves  on  a beach, 
the  waves  become  saturated  at  low  frequencies  first  imposing  a bound 
on  the  peak  energy  density. 

Spectra  of  waves  at  breaking,  or  just  inshore  of  breaking,  are 
shown  in  Figures  10  and  11  plotted  on  a log-log  scale.  The  spectral 
estimates  have  been  normalized  by  dividing  by  the  total  variance  in  an 
effort  to  bring  the  high  end  of  the  spectra  to  a single  line.  The 
spectral  estimates  have  been  block-averaged  over  frequency  bands  such 
that  the  log-energy  values  are  linearly  distributed  on  the  log-frequency 
axis.  The  confidence  intervals  decrease  with  increasing  frequency 
because  of  the  increasing  number  of  frequency  bands  averaged  over.  The 
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confidence  interval  is  noted  on  the  figures.  The  slopes  of  the  log- log 
spectra  at  high  frequencies  for  Wave  Gauges  #1  and  Wave  Gauges  #3  can 
be  approximated  by  a -5  slope.  This  is  indicative  of  a saturation  or 
equilibrium  region  for  wind  generated  waves  in  deep  water.  This  was  a 
surprise,  as  a -3  slope  was  expected  for  the  waves  measured  in  shallow 
water,  Thornton  (1976).  A possible  explanation  is  that  most  of  the  waves 
measured  were  waves  inside  the  initial  break  point  and  would  be  classified 
as  spilling  breakers  or  reformed  waves.  It  may  be  that  these  types  of 
waves  do  not  reach  saturation  condition  as  found  for  the  more  intense 
plunging  and  collapsing  breakers. 

The  normalized  velocity  spectra  are  given  in  Figures  12  and  13.  The 
slopes  of  the  velocity  spectra  at  higher  frequencies  at  Flow  Meter  #1  and 
Flow  Meter  #3  are  both  closely  approximated  by  a -3  slope.  A -3  slope 
is  expected  for  the  equilibrium  region  in  both  deep  and  shallow  water 
and  is  consistent  with  the  wave  slopes. 
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Figure  12.  Onshore  water  particle  velocity  spectra  at 

Flow  Meter  #1 . 


V.  CONCLUSIONS 
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The  values  of  the  pdf's  of  sea  surface  elevation  were  positively 
skewed  verifying  the  asymmetrical  shape  of  the  waves  in  the  breaker  zone. 

On  the  other  hand,  the  skewness  of  the  velocity  components  were  both 
positive  and  negative  and  presented  no  readily  apparent  pattern. 

The  mean  values  of  the  onshore  and  alongshore  components  of  flow  indi- 
cated the  presence  of  a nearshore  circulation  cell  probably  associated  with 
the  rip  current  which  was  frequently  observed  to  the  south  of  the  instru- 
ments. 

The  values  of  the  velocity  energy-density  spectral  components  calculated 
from  wave  spectra  using  linear  theory  indicate  a qualitative,  but  not 
quantitative,  relationship.  In  general,  linear  theory  overestimated  the 
magnitude  of  the  velocity  spectra  because  of  wave  directionality. 

The  high  coherence  between  waves  and  onshore  flow  indicated  that  most 
of  the  motion  in  the  breaking  wave  is  wave-induced.  The  coherence  between 
waves  and  onshore  flow  was  used  to  separate  the  turbulence  to  wave-induced 
velocity  components.  The  measurements  showed  that  over  the  range  of  col- 
lapsing to  spilling  breakers  a reasonable  value  for  the  ratio  of  turbulent 
to  wave-induced  velocity  is  approximately  0.75. 

The  waves  lead  the  onshore  flow  in  spectral  phase  on  the  average 
by  less  than  20  degrees,  implying  that  the  "curling"  crest  of  the  wave 
arrives  prior  to  maximum  water  particle  velocity. 

The  wave  and  velocity  energy-density  spectra  at  higher  frequencies 
exhibited  slopes  of  -5  and  -3,  respectively.  This  is  indicative  of  deep 
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water  wave  conditions.  A possible  explanation  for  the  surprising  spectral 
slope  is  that  most  of  the  waves  measured  were  either  spilling  breakers  or 


reformed  waves, 
tion  conditions 


It  may  be  that  the  measured  waves  did  not  attain  satura- 
during  breaking,  as  do  plunging  and  collapsing  breakers. 


APPENDIX  B 

CAPACITANCE  WAVE  GAUGE  CALIBRATIONS 


DEPTH  OF  miERSION  (FEET) 


WG  #3 
FM  #3 
ONSHORE 


ALONGSHORE 


♦Calibration  Additive  Factor  for  wave  heights  is  given  in  meters; 
GALA  for  flow  is  given  in  meters  per  second. 

♦♦Calibration  Multiplication  Factor  for  wave  heights  is  given  in 
meters  per  bit;  CALM  for  flows  is  given  in  meters  per  second  per  bit. 
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APPENDIX  E 

POWER  COHERENCE  AND  PHASE  SPECTRA 

WAVE  GAUGE  #1  - Onshore  Flow  at 
FLOW  METER  #1 
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